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I. INTRODUCTION 

The correct understanding of the energy, momentum, 
and angular momentum in gravity theories is a prominent 
physical problem that has a long and complicated his¬ 
tory. In order to solve this problem, numerous formalisms 
were developed to derive what are generally known as 
the conservation laws. There are two classes of conser¬ 
vation laws in any gravity theory. One class of conser¬ 
vation laws is formulated solely in terms of the dynam¬ 
ical variables that describe the gravitational field itself. 
These variables characterize the geometry of spacetime 
without involving additional structures of physical (non- 
geometrical) nature. The resulting conserved charges do 
not have tensor transformation properties under general 
coordinate transformations. 

There is, however, another class of conservation laws 
that lead to conserved charges and that turn out to be 
true scalars. In deriving such conservation laws, one 
usually has to deal with, besides the gravitational field 
variables, additional physical structures such as vector 
fields. As it is well known, vector fields generate dif- 
feomorphisms on a spacetime manifold. What is more 
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important, one can associate with a vector field Cf a cur¬ 
rent that is conserved under some conditions. In Ein¬ 
stein’s general relativity theory this can be illustrated as 
follows. 

Given a covariantly conserved symmetric energy- 
momentum tensor, ViTjU = 0, we find (contracting with 
C,^) a relation: V^J* * * § = Here we defined the 

current J* := furthermore, Vj is the Riemannian 

covariant derivative and Cq is the Lie derivative along 
the vector field . If is a Killing vector, the current 
is conserved = 0. This fact establishes a remark¬ 
able relation between the symmetries of the spacetime 
and the conserved currents generated by these symme¬ 
tries. Physically, the vector field is usually related to the 
reference frame motion of an observer. 

In this paper we extend this observation to a general 
framework of gravity theories that encompass possible 
non-Riemannian geometries and nonminimal couplings 
of matter to the gravitational field. Theories of this kind 
attracted considerable attention in the literature. 

The history of the construction of conserved quanti¬ 
ties in relativistic gravity is rich and long. The earliest 
relevant construction is the Komar charge and its close 
relatives [1-3]. Essential contributions to this approach 
include [4-19]. Most of these works were confined to the 
purely Riemannian geometrical framework of Einstein’s 
general relativity. Extensions to more general geometries 
were studied in [20-26] and more recently by us [27-29] 
and [30-32]. 
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The structure of the paper is as follows. In the next 
Sec. II we give an overview of the basic geometrical no¬ 
tions and operations which underlie our study. Sec. Ill 
presents the essentials of generalized gravity theory. We 
give a detailed derivation of the currents associated with 
the vector fields (diffeomorphisms) on the spacetime 
manifold, and specify the conditions under which these 
currents are conserved. In particular, we demonstrate the 
importance of generalized Killing vector fields, the prop¬ 
erties of which we discuss in Sec. IV. In Sec. V we use 
our findings to obtain the generalized conserved current 
in metric-affine gravity (MAG) with a possible nonmin- 
imal coupling. Finally, we also show that the equations 
of motion for extended microstructured test bodies also 
admit a conserved quantity that is closely related (and 
has a similar structure) to the conserved currents derived 
in the field-theoretic framework. We conclude our paper 
in Sec. VI with a discussion of the results obtained and 
with an outlook of their possible applications. 

Our notations and conventions are summarized at the 
end of the paper, in Appendix A. 

II. PRELIMINARIES: SPACETIME GEOMETRY 

In the metric-af&ne theory of gravity, the gravitational 
physics is described by the dynamics of the geometri¬ 
cal structure of spacetime. The latter is encoded in two 
fields: the metric tensor gij and an independent linear 
connection T . The latter is not necessarily symmetric 
and/or compatible with the metric. From the geometrical 
point of view, the metric introduces lengths and angles of 
vectors, and thereby determines the distances (intervals) 
between points on the spacetime manifold. The connec¬ 
tion introduces the notion of parallel transport and de¬ 
fines the covariant differentiation Vfe of tensor fields. 

Locally, a spacetime diffeomorphism can be described 
as a small translation in the spacetime manifold, which 
technically is represented by the variation of the space- 
time coordinates 

Sx^ = ee*(a;). (1) 

Here e is an infinitesimal constant parameter and e®(a;) 
is an arbitrary, but finite, vector field. Under the action 
of diffeomorphisms, the geometrical variables transform 
as 

Sgzj = -e{d,£^)gkj - e{dj£^) g,k, ( 2 ) 

STk^ = -e{dk£^)Ti^ - eid,£^)ru^ 

+ e{di£^)rk/ - edl,£^. (3) 

In general, the geometry of a metric-affine manifold is 
exhaustively characterized by three tensors: the curva¬ 
ture, the torsion and the nonmetricity. They are defined 
[33] as follows: 

Rki^ ■■= dkTi^ - diTk^^ + - TiJTk,’-, ( 4 ) 

Ffei* - (5) 

Qkij •— ^k9ij — T ^ki 9lj ^kj Qil- (6) 


The curvature and the torsion tensors determine the 
commutator of the covariant derivatives. For a tensor 
of arbitrary rank and index structure: 

(VaVf,- Vf,Va)A"i -"'=di...d, = 

k 

I 

i=i 

The Ricci tensor is introduced by Rij := Rkij^, and the 
curvature scalar is i? := g'^^Rij. 

A general metric-affine spacetime {Rkii^ 7 ^ 0, Tu 7 ^ 0, 
Qkij 7 ^ 0 ) incorporates several other spacetimes as spe¬ 
cial cases. The Riemannian connection T^j^ is uniquely 
determined by the conditions of vanishing torsion and 
nonmetricity which yield explicitly 

Tfc/ = ]^g'’\djgki + dkgij - digkj). (8) 

Here and in the following, a tilde over a symbol denotes 
a Riemannian object (such as the curvature tensor) or a 
Riemannian operator (such as the covariant derivative) 
constructed from the Christoffel symbols (8). The devi¬ 
ation of the geometry from the Riemannian one is then 
conveniently described by the distortion tensor 

iVfc/ := Ffc/ - Ffe/. (9) 

The definitions (5) and (6) allows to find the distortion 
tensor in terms of the torsion and nonmetricity. Explic¬ 
itly, 

= — —(Tfcj* -|- T^kj + T^jk) + — Qkj^ — Qjk)- 

( 10 ) 

Conversely, one can use this to express the torsion and 
nonmetricity tensors in terms of the distortion, 

Tfe/ = -2iV[fe,f, (11) 

Qkij = —‘2Nk{i’‘gj)i- ( 12 ) 

Substituting (9) into (4), we find the relation between the 
non-Riemannian and the Riemannian curvature tensors, 

RadJ’ = Radc^-^aNdc’^+VdNaJ’+Nan’^Ndc^-Ndn^Nac^ 

(13) 

Applying the covariant derivative to (4)-(6) and antisym¬ 
metrizing, we derive the Bianchi identities [33]: 

= R[klnf + T[kl'^Tn]rn 
^ [nQk]ij ~ Rnk{ij) ■ 


(14) 

(15) 

(16) 
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A. Tensors, densities, and covariant differential 
operators 

Along with tensors, an important role in physics is 
played by densities. A fundamental density ^/—g is con¬ 
structed from the determinant of the metric, g =detgij. 
Under diffeomorphisms (1) it transforms as 

Sy/^ = (17) 

This is a direct consequence of (2). From any ten¬ 
sor - one can construct a density ’ = 

yJ—gBi,,,^ '. Although in this paper we will encounter 
only such objects, it is worthwhile to notice that not all 
densities are of this type, since they can have different 
weights. The fundamental density \/—g and all other 
densities discussed here have weight -|-1. See the exhaus¬ 
tive presentation in the book of Synge and Schild [34]. 

There are two kinds of covariant differential operators 
on the spacetime manifold, depending on whether the 
connection is involved or not. The Lie derivative Cg is 
defined along any arbitrary vector field e® and it maps 
tensors (densities) into tensors (densities) of the same 
rank. Let us recall the explicit form of the Lie derivative 
of the metric and the distortion: 

Cegij = e’^dkgij + {d^e'^)gkj + {dje'")gik, (18) 

ndje^Nkn - (9ne*)iVfc/. (19) 

In contrast, a covariant derivative Vfe raises the rank 
of tensors (densities) and it is determined by the linear 
connection Tfcj®. Moreover, there are different covariant 
derivatives which arise for different connections that may 
coexist on the same manifold. 

A mathematical fact is helpful in this respect: every 
third rank tensor Xkj^ defines a map of one connection 
into a different new connection 

Tfe/ ^ Tkj^+Xk/. (20) 

There are important special cases of such a map. One 
example is obtained for Xi~/ = Nk/: then the connection 
Tfcj® is mapped into the Riemannian Christoffel symbols, 
Lfcj® = Tfcj® -I- Nk/, in accordance with (9). 

Another interesting case arises for Xkj^ = Tjk^. The 
result of the mapping 

Ffc/ = Tfe/ + TjU = Tfe/ + TjU - Tfc/ = r,U (21) 

is then called a transposed connection^ or associated con¬ 
nection, see [35, 36]. 

The importance of the transposed connection is man¬ 
ifest in the following observation. Although the Lie 
derivative is a covariant operator - this is not apparent 
since it is based on partial derivatives - one can make 
everything explicitly covariant by noticing that it is pos¬ 
sible to recast (18) and (19) into equivalent forms 

^egij = £^Xkgij + {X ie^)gkj + {X j£^)gik, ( 22 ) 

C,Nki = e^XnNk/ + (Vfce”)iV„/ 

+ {Xj£^)Nkn^ - (V„e®)iVfe/. (23) 


By the same token we can “covariantize” the Lie deriva¬ 
tives for all other tensors of any structure and of arbitrary 
rank. 

A more nontrivial (and less known) fact is that we can 
define the Lie derivatives also for objects which are not 
tensors. In particular, the Lie derivative of the connec¬ 
tion then reads [35]: 

C.Tki = e^diVk/ + {dk£^) T,/ + {d.e^) Tfe,® 

- idie^)Tk/+ (24) 

^XkXjE^ - RkiieK (25) 

This quantity measures the noncommutativity of the Lie 
derivative with the covariant derivative 

ic,Xk-XkCM"^-"’‘d^...d, 

k 

I 

-Y,i^,TkdM^^-^'‘d....b...dr (26) 

J-1 

The connection Tfcj®, the transposed connection T^j®, 
and the Riemannian connection Tfcj® define the three re¬ 
spective covariant derivatives: Vfc, Xk, and Xk- 

We will assume that these differential operators act on 
tensors. In addition, we will need the covariant operators 
that act on densities. For an arbitrary tensor density 
18®®i...'^ ’ we introduce the covariant divergence 

v„$®®,„/ - ■■■■ 

(27) 

which produces again a tensor density. We denote a simi¬ 
lar differential operation constructed with the help of the 
Riemannian connection by 

- := - - 

(28) 

When -, we find 

- = v^V^B®®,../-, (29) 

V„S®®,.U■■■■ = (30) 

where we introduced a modified covariant derivative 

V*:=V, + Ah'=. (31) 

A remark explaining our notation is in order. As one 

knows, a tensor is special case of a density with zero 
weight. So, strictly speaking, the introduction of the 
numerous superscript accents above may be viewed as 
something redundant. However, we find it more conve¬ 
nient to explicitly distinguish densities from tensors with 
the help of the “Fraktur” font for the former, and the 
“Roman” font for the latter. Accordingly we use differ¬ 
ent accented symbols for the derivatives. In particular, 
one should note that (31) acts on tensors, in contrast to 
(27) and (28) which act on densities. 
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B. Matter variables 

We will not specialize the discussion of matter to any 
particular physical field. It will be more convenient to 
describe matter by a generalized field The range 
of the indices A, B,... is not important in our study. 
However, we do need to know the behavior of the matter 
field under spacetime diffeomorphisms. We assume that 
under the transformation (1), these fields satisfy 

= - €{die^) ■ (32) 

Here {(t"^b)/ are the generators of general coordinate 
transformations that satisfy the commutation relations 

= ia^ByS^-{a^B),'^6l (33) 

We immediately recognize in (33) the Lie algebra of the 
general linear group GL(A,R). This fact is closely re¬ 
lated to the standard gauge-theoretic interpretation [37] 
of metric-affine gravity as the gauge theory of the general 
affine group G'H(4, i?), which is a semidirect product of 
spacetime translation group times GL{A,R). 

The transformation properties (32) determine the form 
of the covariant and the Lie derivative of a matter field: 



:= dki>^ - Tki^ {a^B)/i’^, 

(34) 


:= + {dy){a^B)/ 

(35) 


= e'=VfcV'^ + (V.e^)((7^i3)/^^. 

(36) 


The commutators of these differential operators read 

(VfcVz - v,Vfe)V'^ = - Rkyy^Bhy^ 

(37) 

(/:,Vfc - = - {c,Ty){a^By^^- (38) 


III. METRIC-AFFINE GRAVITY: FIELD 
EQUATIONS AND CURRENTS 


renorinalizable theory of the quantized gravitational field 
[44-48]. 

In the standard formulation of MAG as a gauge the¬ 
ory [37], the gravitational gauge potentials are identified 
with the metric, coframe, and the linear connection. The 
corresponding gravitational field strengths are then the 
nonmetricity, the torsion, and the curvature, respectively. 

It will be convenient to describe all the dynamical vari¬ 
ables - including the gravitational (geometrical) and ma¬ 
terial fields - collectively by means of a multiplet, which 
we denote by = (g^, V'"^)- The range of the 

multi-index J is not important at present and will be 
specified when needed. 

A. Lagrange-Noether analysis 

The dynamics of the interacting gravitational and mat¬ 
ter fields is determined by a general action 

(39) 

The Lagrangian density £ depends arbitrarily on its ar¬ 
guments. Let us investigate the variation of the action 
under transformations of the spacetime coordinates and 
the fields which, quite generally, read as follows: 

a;'® (a;) = + 5a;®, (40) 

^'yx') = ^yx)+S<S>yx). (41) 

Then, using the substantial variation defined by 
5(s)<I’'^ := $''^(a:) — 4’'^(a:) = 5$'^ — Sx^dk^"^, we derive 
in a standard way the total variation of the action 

5/ = y yx [5(,)£ -b (£5a:®)] . (42) 

Assuming the invariance of the action under (40) and 
(41), we find the so-called Lie differential equation 

5(,)£-b9* (£5a;®) = 0. (43) 


The explicit form of the dynamical equations of the 
gravitational field is irrelevant for the conservation laws 
that will form the basis for the derivation of the test body 
equations of motion. However, for completeness, we dis¬ 
cuss here the field equations of a general metric-affine 
theory of gravity. The standard understanding of MAG 
is its interpretation as a gauge theory based on the gen¬ 
eral affine group GA(4, R), which is a semidirect product 
of the general linear group GL(A, R), and the group of lo¬ 
cal translations [37]. The corresponding gauge-theoretic 
formalism generalizes the approach of Sciama and Kibble 
[38, 39]; for more details about gauge gravity theories, see 
[40-43]. Besides its many interesting properties, MAG of¬ 
fers the possibility of a unification of gravity with other 
physical interactions on the same gauge-theoretic prin¬ 
ciples, and contributes to the solution of the quantum 
gravity quest with encouraging attempts to construct a 


With the help of the chain rule we can write 




d£ 

d{d^ 




(44) 


and using the commutativity of the substantial variation 
with the partial derivative, S(^s)di = 9i5(s), we recast (43) 
into a balance equation 


^ 3 ,., 4 ^ + 8 .(£&- + 


d£, 






0. (45) 


Here we denote the variational derivative 

5 £ 5 £ ^ f d£ 

S<^J ■“ d<i>J ~ " \d{di<^y 

as usual. When the variational derivatives are put equal 
to zero, we obtain the system of classical field equations. 
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The configurations of variables that satisfy 6£,/6^'^ = 0 
are called “on-shell”. 

The balance equation (45) is an identity, i.e. it is valid 
for all configurations of the gravitational and matter 
fields irrespectively of the fact that they are “on-shell” or 
“off-shell”. Equation (45) gives rise to various identities 
and conservation laws which we derive in this paper. 

We can apply (45) to different symmetries of the physi¬ 
cal system under consideration. Of particular importance 
is the diffeomorphism (or general coordinate) invariance 
of the action (39). Substituting (1), (2), (3), and (32) 
into (45), we recast the latter (after dropping the overall 
infinitesimal constant e) into 

Qk d^e’^ + iofcb 52 ^ gyk ^ (47) 

The functions (with n = 0,1,2,3) are deter¬ 

mined by the Lagrangian of the theory. Their explicit 
form is given in Appendix B. 

In view of the arbitrariness of the vector field and 
their derivatives, we obtain a set of Noether identities 

= 0 (48) 


for n = 0,1, 2, 3. 


B. Generalized current 

Every vector field on the spacetime manifold induces 
a current which is conserved under certain conditions. 
We can derive it with the help of the balance equation 
(45) as follows. Let us consider a map of the manifold 
(diffeomorphism) induced by a vector field, as in Eq. (1). 
Locally, such a diffeomorphism acts on the gravitational 
and matter variables = {gij ,Tki^ ,ip^) by means of 
the Lie derivatives so that 


= -eC.Tkj^ 


(49) 

(50) 

(51) 


Inserting this into (45) and dividing by the infinitesimal 
parameter e, we can define the current density 


d := 


dd^gki 

d£ 


^e9kl + 




dd^Tkn^ 




e-L kn 


dditp^ 

and observe that it satisfies 

Old — e ^e9kl 

09 kl 0l kn 




(52) 


-UT 


kn 




Sip 


(53) 


Using (18), (25), and (36), we can display the structure 
of the current as follows 


a* = + (VfeV^e™)^, 


kni 


(54) 


Here we introduced the densities 


dm := 


dZ 




ddiip^ 

^mkl cvrx yn „ ^kml i 


^^iQkl 

■■= - 9lr. 


ddZT 

dZ 


kl 


u^lB 


kni 


dZ 


ddi'ip^ 


{(T B)m^1p 


kn 


(55) 

(56) 

(57) 


C. Current and superpotential 


Expanding the Lie derivatives in (52), we can recast 
the current into a different equivalent form. Namely, a 
lengthy direct computation yields 




— e-' 


'5Z 5Z , A ^ 


- (ymen 


L'-L rrr). 


(58) 


Here we introduced a new density 


0 := - e‘ 


- 9lk 


dZ 


ij.B 


B)k'4) 


ddjgu'’"' ' ddjtp^ 

5£ . a£ a£ 

-L kn “r rx o t-, ^ ^ nk 


dTji'^ ddjTin 

dZ , dZ 

^ nm " 


hd T ^ 


f)f) V ^ 

^ nm 


- [dnS'") 


dZ 


J 

dZ 


ddjVin^ dd(rPZj)i'' 


(59) 


The first line on the right-hand side of (58) vanishes for 
the diffeomorphism invariant (generally covariant) theo¬ 
ries in view of the Noether identities (48). As a result, we 
obtain “on-shell” (i.e., when 6Z/dgij = 0, 6Z/STkj^ = 0, 
and SZ/Stp"^ = 0): 


(60) 

This means that on the classical field equations, the cur¬ 
rent d® is derived from the superpotential density . 


D. Explicitly covariant current 

It is not obvious that the complicated expressions (52), 
(58), and (59) are truly covariant objects under the ac¬ 
tion of diffeomorphisms. A direct demonstration is pos¬ 
sible, but it is somewhat long to present it here. Instead, 
we will show the covariance by specifying the form of 
the Lagrangian. Intuitively it is clear (and can be rig¬ 
orously proven) that the Lagrangian of a generally co¬ 
variant theory should be a function of covariant objects. 
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This means, in particular, that the derivatives of the ba¬ 
sic variables , 1 !;^) should only appear in 

the form of explicitly covariant combinations. In simple 
terms. 


We then immediately compute the partial derivatives 
d£ dZ 


ddjgu 

d£ 

ddjVin'^ 

d£ 


= 2 


dQjii ’ 

d£ 


dRj, 

d£ 


d£ ^ d£ 

— 2 „ „ gji 


d£ 




dQkii' 
d£ 


+2 


dRkmJ 


-r„ 


dZ 


^dTk^i dVki’^ 

I g r 


OR: 


mki 


(62) 

(63) 

(64) 

(o-'^B)jV^ 

(65) 


mj 


Substituting these expressions, we obtain the explicitly 
covariant current 


d* = -e^ 


'5Z ^Z , A . ^ ^ 


- (V^£") 


(5£ 

SVim' 




( 66 ) 


and the explicitly covariant superpotential 

_ op 

i?-=2e'=^+2(V„e'=)^. (67) 

L/ -t 7 7 UlL-j/ir), 




In our work [27], we studied the definition of conserved 
currents within models with local Lorentz invariance and 
vanishing nonmetricity. In those models, one can define 
different conserved currents depending on the choice of 
a so-called generalized Lie derivative. See [27, 29] for 
further details. On the other hand, the current (66) 
and the superpotential (67) are defined for any metric- 
affine model. Comparing both results for the restricted 
case of vanishing nonmetricity, one can verify after some 
straightforward algebra, that the superpotential (67) re¬ 
duces to the components of the potential 2-form H of 
[27], provided one uses the Yano choice for the general¬ 
ized Lie derivative [49—51]. Therefore, the corresponding 
currents coincide when the classical field equations for 
the total system are satisfied. As a consequence, in the 
case of the Hilbert-Einstein Lagrangian, the conserved 
quantities reduce, in vacuum, to the well known Komar 
charges [1, 2]. 


E. Conservation of the general current 

The current 3* defined in (52), or equivalently in (66), 
is conserved when the right-hand side of the balance 
equation (53) vanishes. 


This is the case, for instance, for the “on-shell” field con¬ 
figuration when both the gravitational and matter vari¬ 
ables satisfy the classical field equations 6Z/6gij = 0, 
SZ/STkj'^ = 0, and SZ/Stl;^ = 0. Alternatively, the cur¬ 
rent is conserved, even “off-shell”, when the Lie deriva¬ 
tives of all field variables vanish along a particular choice 
of the vector field £*. In these cases, when the vector field 
is not arbitrary, but some particular special case, we de¬ 
note it by C, so that it satisfies C(^gij = 0, C<;Zkj'' = 0, 
and = 0. 

The last assumptions are, however, a bit too strong. 
A milder condition for the conservation of the current is 
the “on-shell” ansatz for the matter fields, dZ/Stp^ = 0, 
combined with the vanishing of only the Lie derivatives 
for the gravitational fields along a vector field Q-. Cc^gij = 
0 and C(iZkj^ = 0. We discuss the geometrical meaning 
of the latter conditions in Sec. IV. 

Note also that (68) can be written in terms of the co¬ 
variant derivatives V and V since, by the general defini¬ 
tions (27) and (28), we have 


ya* = (69) 

IV. SYMMETRIES IN MAG: GENERALIZED 
KILLING VECTORS 

As is well known, symmetries of a Riemannian space- 
time are generated by Killing vector fields. Each such 
field defines a so-called motion of the spacetime mani¬ 
fold, that is a diffeomorphism which preserves the metric 
dij ■ 

Suppose C* is a Killing vector field. By definition, it 
satisfies 


V,0+V,C. = 0. (70) 

By differentiation, we derive from this the second covari¬ 
ant derivative 


]Ck = Rjki Cl- ( 71 ) 

We apply another covariant derivative and antisym¬ 
metrize: 


V[„Vj]VjCfc = y [n{R\jk\i\\i) ■ (72) 

After some algebra, the last equation is recast into 

C^VnRijkl + Rnjkl^ + Rinkl^ jC^ 

+R^JnlVkC + R^kn^lC = 0. (73) 

Equations (70), (71), and (73) have a geometrical 
meaning: 


II 

0 

(74) 

= 0, 

(75) 

^cRijkl = 0. 

(76) 


= 0 . 


( 68 ) 
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That is, the Lie derivatives along the Killing vector field ( 
vanish for all Riemannian geometrical objects. Moreover, 
one can show that the same is true for all higher covariant 
derivatives of the Riemannian curvature tensor [49]: 

Cc{Vn,...Vn^R,jkl) =0- (77) 

It is worthwhile to mention that in the Riemannian 
framework of Einstein’s general relativity, one can define 
various symmetries generated by the vector fields (dif- 
feomorphisms). For example, if ( does not satisfy (74) 
but fulfills (75), such a vector field is not a Killing but a 
so-called called affine collineation. Alternatively, if both 

(74) and (75) are not true but C is characterized by the 
property (76), it is called a curvature collineation. Along 
with the standard Killing vector fields, these new fields 
contain important information about the symmetries of 
the spacetime, and it is possible to use them to define ad¬ 
ditional conservation laws. A comprehensive discussion 
of such symmetries (and related ones such as homoth- 
etic, conformal, projective, Ricci collineations etc) can 
be found in [52, 53], and different applications are stud¬ 
ied in [54-56], for example. 

Let us generalize the notion of a symmetry to the 
metric-afhne spacetime. We take an ordinary Killing vec¬ 
tor field C and postulate the vanishing of the Lie deriva¬ 
tive 

CcNk/ = 0 (78) 

of the distortion tensor. Combining this with (10) and 

(75) , we find an equivalent formulation 

Ccg^j = 0, (79) 

= 0. (80) 

We call a vector field that satisfies (79) and (80) a gen¬ 
eralized Killing vector of the metric-affine spacetime. By 
definition, such a C generates a diffeomorphism of the 
spacetime manifold that is simultaneously an isometry 
(79) and an isoparallelism (80). 

Since the Lie derivative along a Killing vector com¬ 
mutes with the covariant derivative, LqS/ i = 'S /see 
(26), we conclude from (13) and (76) that the general¬ 
ized Killing vector leaves the non-Riemannian curvature 
tensor invariant 

C^Rkij^ = 0. (81) 

It is also straightforward to verify that 

(V„i ... = 0 (82) 

for any number of covariant derivatives of the curvature. 

Finally, combining (11) and (12) with (74) and (78), 
we verify that 

= 0, (83) 

C-QQkij = 0. (84) 


Summarizing, the Lie derivative of all main geometrical 
objects vanishes along a generalized Killing vector field. 

Later we will show that generalized Killing vectors 
have an important property: they induce conserved 
quantities on the metric-affine spacetime. 

V. MODIFIED MAG MODELS WITH A 
POSSIBLE NON-MINIMAL COUPLING 

A. Gravitational field dynamics 

The general formalism which we developed in Sec. Ill 
will now be applied to specific models of the gravita¬ 
tional field coupled with matter. However, before we 
go into technical details, the following remark is in or¬ 
der. Strictly speaking, the metric-affine framework is ap¬ 
plicable not only to generalized gravity theories which 
are based on the non-Riemannian spacetime geometries, 
but it is equally useful for the study of Einstein’s GR 
and its modifications in a purely Riemannian geometri¬ 
cal context. Technically, this requires the introduction 
of additional variables which play the role of Lagrange 
multipliers, and this changes the physical meaning of the 
sources in the field equations. One should take note of 
this subtlety. 

Here we use the formulation of MAG, in which grav¬ 
ity is described by the set of fundamental field variables 
which consists of the independent metric gij and connec¬ 
tion Tfci-l. Such an approach was developed in [57-66], 
and this is alternative to the formalism which includes 
also a coframe field [37] as a gravitational field variable. 
It is instructive to compare the field equations in the 
different formulations of MAG, and in particular, it is 
necessary to clarify the role and place of the canonical 
energy-momentum tensor as a source of the gravitational 
field. Since one does not have the coframe (tetrad) among 
the fundamental variables, the corresponding field equa¬ 
tion is absent. Here we demonstrate that one can always 
rearrange the field equations of MAG in such a way that 
the canonical energy-momentum tensor is recovered as 
one of the sources of the gravitational field. 

For a large class of MAG models, the total Lagrangian 
of interacting gravitational and matter fields reads 

£ = + (85) 

In general, the gravitational Lagrangian is constructed as 
a diffeomorphism covariant density function of the cur¬ 
vature, torsion, and nonmetricity, 

03 = QJ [g^j , \Tki\Qkii), (86) 

whereas the matter Lagrangian depends on the matter 
fields and their covariant derivatives (34): 

£mat = £mat(5*J,V^^,V,V'^). (87) 

In the current literature a considerable attention is 
paid to the study of so-called modified models in which 


one assumes a possibility that the gravitational field may 
interact with matter in a nonminimal way. Accordingly, 
we here also allow for nonminimal interaction of the mat¬ 
ter to the gravitational field via the coupling function 

F = F{gij,Rijk\Tki^ ,Qkij). ( 88 ) 


On the right-hand side, the matter source is now repre¬ 
sented by the metrical energy-momentum tensor which 
is defined by 

ftp 

b,:=2^. (100) 


When = 1, we recover the minimal coupling case. 
Let us write down the field equations of metric-affine 
gravity for that case. This can be done in several equiv¬ 
alent ways. The standard form is the set of the so-called 
“first” and “second” field equations. Using the covariant 
derivative for densities defined by (27), the field equa¬ 
tions are given by 

= -TU, (89) 

= ©*/. (90) 

Here the generalized gravitational field momenta densi¬ 
ties are introduced by 






dRkii^ ’ 




dQ kij 


(91) 

(92) 

(93) 


and the gravitational hypermomentum and the general¬ 
ized energy-momentum densities are constructed as 


•yki 


3^ 




+ Tkl^Sj^‘n+Fkln”^:fj 


m cLiln 


(94) 

(95) 


The sources of the gravitational field are the canonical 
energy-momentum tensor density and the canonical hy¬ 
permomentum density of matter, respectively: 


TU := 


:= 


d&n 


9Umat 

^^k^^ 


9£mat 




(cr B)jV 


B 


(96) 

(97) 


The usual spin density arises as the antisymmetric part 
of the hypermomentum. 


T' • © 


ij] 


(98) 


whereas the trace ©^ = ©*i^ is the dilation current den¬ 
sity. The symmetric traceless part describes the proper 
hypermomentum [37]. 

It is straightforward to verify that instead of the first 
field equation (89), one can use the so-called zeroth field 
equation which reads 


2—= LU (99) 

^9ij 


The system (89) and (90) is completely equivalent to the 
system (99) and (90), and it is a matter of convenience 
which one to solve. 


B. Conservation laws 

The dynamics of extended bodies in metric-affine 
spaces can be derived by integrating the conservation 
laws. The latter are obtained from the Noether iden¬ 
tities. After some algebra, the identities (48) are recast 
into the following set of conservation laws 

FTfc* = FtU + V„ (U©\,"), (101) 

V, (UTU) = F - &^n'Rkim^ - le^Qk^J^ 

— UmatVfeF. (102) 

These relations hold “on-shell” when the matter variables 
satisfy the classical field equations, and are valid for 
general nonminimal coupling. 


C. Rewriting the conservation laws 

Using (31) and decomposing the connection into the 
Riemannian and non-Riemannian parts, c.f. Eq. (9), we 
can recast the conservation law (101) into an equivalent 
form: 

^jiF&k^ = Fflk^ - tfc* + - Nr.k'^&m^. 

(103) 

In a similar way we can rewrite the conservation law 
(102). At first, with the help of (11) and (12) we no¬ 
tice that 

F - ^e^Qk^J^ = F{if - %^)Nk/ + F%^N,k'. 

(104) 

Then substituting here (103) and making use of (31), (9), 
and the curvature decomposition (13), after some algebra 
we recast (102) into 

V,{F(Tfe^ + e^n^Nkm^} 

= -U©™„*(Rfe,™"-VfciV,™") 

- UmatVfeF. (105) 

For the minimal coupling case, such a conservation law 
was derived in [27, 67]. The importance of this form of 
the energy-momentum conservation law lies in the clear 
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separation of the Riemannian and non-Riemannian geo¬ 
metrical variables. As we see, the post-Riemannian ge¬ 
ometry enters (105) only in the form of the distortion ten¬ 
sor which is coupled only to the hypermomentum 
current density ©"*„*. This means that, in the minimal 
coupling case, ordinary matter - i.e. without microstruc¬ 
ture, ©'"n* = 0 - does not couple to the non-Riemannian 
geometry. In contrast, in the nonminimal coupling case, 
the derivative of the coupling function F on the right- 
hand side of (105) may lead to a coupling between non- 
Riemannian structures and ordinary matter. 


D. Conserved current induced by a spacetime 
symmetry 


As we have shown, every vector field generates a con¬ 
served current. It is straightforward to show that for a 
Lagrangian density of the form £ = A'£mat j as defined in 
(87) and ( 88 ), if we consider a generalized Killing vector 
field the current (54) reads, 

y := F - (V^r )e™n*] . (106) 

Here we have used the definitions (96) and (97), and the 
conditions (79) and (80). 

It is instructive to derive this result directly from the 
field equations. Namely, let us contract equation (103) 
with and contract equation (105) with and then 
subtract the resulting expressions. Note that the contrac¬ 
tion = 0 vanishes because the first factor is a 

symmetric tensor and the second one is skew-symmetric. 
Then after some algebra we find 

= Fe^^n^CcN^m^ - H^^tCcF (107) 


The right-hand side of (107) depends linearly on the Lie 
derivatives along the Killing vector: and 

see (19) and (23). 

When is a generalized Killing vector, we have 
= 0 in view of (78). Furthermore, recalling that 
F = F{gij,RMj\Tki^,Qkij), we hnd 


dF OF 

OF dF 

= 0 , (108) 


dT, 


kj 


dQ kij 


by making use of (74), (81), (83), and (84). 

As a result, the right-hand side of (107) vanishes for 
the generalized Killing vector field, and we conclude that 
the induced current (106) is conserved. 

This generalizes the earlier results reported in [27, 28, 
36]. In Sec. VF we will show that there is a conserved 
quantity constructed from the multipole moments which 
is a direct counterpart of the induced current (106). It is 
worthwhile to give an equivalent form of the latter: 


y = F 


+ Q'^n^Nkm^) - {ymC)& 


(109) 


E. Test body equations of motion 

The details of the multipole approximation scheme are 
given elsewhere; see Appendix C for the definitions of the 
integrated moments. Here we limit ourselves to the pole- 
dipole equations of motion. Following [ 68 ], we introduce 
the total orbital and the total spin angular moments 

L'^b 2p[ab] ^ _ _2/j[a6], 

where p°‘^ and h°‘^ are generalized multipolar moments. 
The generalized total energy-momentum 4-vector and the 
generalized total angular momentum are given by 

:= T(p“ + N‘^cdh^‘^) + p’’‘^VbF, (111) 

■= F{L^^ + S^^). (112) 

Then the pole-dipole equations of motion take the form 

1 ~ ~ 

+ Fq^’’‘^V‘^Ndab 

as 2 

- (113) 

D 'T°‘b 
as 

_ 2 ^[ayb]^^ (114) 

Here is the 4-vector of velocity of the test particle, 

= u“Vq, with s the proper time of the particle, 
whereas < 7 “*''^, ^ and are integrated moments relevant 
for the description. See App. C, as well as [ 68 ] for further 
details. 


F. Conserved quantity for extended test bodies 

The equations of motion for the multipole moments 
are derived from the conservation laws of the energy- 
momentum and the hypermomentum currents and 
©"*„*. In Sec. VD we have demonstrated that every 
generalized Killing vector induces a conserved current 
depending on these two quantities, see Eqs. (106) and 
(109). Quite remarkably, there is a direct counterpart 
of such an induced current built from the multipole mo¬ 
ments. 

Let be a generalized Killing vector, and let us con¬ 
tract equation (113) with C,a and equation (114) with 
and then take the sum. This yields 

£ (^“Ca + = Fq^^^jr^Ndcb 

- ^C^F-C^XcF.ill5) 

On the right-hand side of (115) the Lie derivatives of the 
distortion tensor and of the coupling function vanish in 
view of (78) and (108). 

Consequently, we conclude that for every generalized 
Killing vector field the quantity 

7’“Ca + = const. 


(116) 
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is conserved along a trajectory of an extended body. 

We thus observe a complete consistency between (106) 
and (116), as well as between (107) and (115). 

The conserved quantities with the similar structure in 
(116) were derived previously for extended test bodies 
with mass and spin in Einstein’s general relativity [69- 
71] and in Einstein-Cartan gravity [36]. Their various 
applications were recently discussed in [72, 73]. 


G. Relation between theories of fields and particles 


The similarity between the currents in a field-theoretic 
picture and the integrals of motion in particle mechanics 
is not occasional. Here we demonstrate that the current 
(106), (109) actually generates the conserved quantity 
(116) in the dipole approximation. Moreover, we will 
show that there exists a generalization of (116) to an 
arbitrary multipole order. 

As it is well established, the equations of motion of 
extended test bodies in external gravitational and other 
classical fields are derived from the corresponding con¬ 
servations laws, see [74-77] for example. In simple 
terms, the dynamical description of an extended body 
is achieved by assigning to a body a set (infinite, in gen¬ 
eral) of multipole moments. The latter are determined 
by certain integrals of the conserved currents (energy- 
momentum tensor and other Noether currents) over the 
body. The multipole moments represent such character¬ 
istics of the body as its total mass, charge, momentum 
etc. By integrating the partial differential conservation 
laws, one obtains a system of ordinary differential equa¬ 
tions of motion for the moments. 

Here we apply the covariant multipole expansion of [76] 
to the current (106). The method is based on Synge’s [78] 
notion of the “world function” cr, which introduces a co¬ 
variant generalization of the finite distance between the 
spacetime points x and y. Basic definitions and notation 
are summarized in Appendix C. By construction, this ob¬ 
ject (t(x, y) depends on two arguments, and as a result 
the Synge formalism deals with arbitrary bitensor densi¬ 
ties The most important technical 

tool is represented by the lemma [76]: 


£ J 

Y.(s) 


E(s) 




E(s) 


We introduce the 4-velocity := dxy^/ds, with the 
proper time s, and denote ^ the integrals are 

performed with an arbitrary bitensor density (a;, y) 

over an arbitrary spatial hypersurface E. See [76] for 
more details on the integrals and the construction of . 
As in all our previous papers, we use a condensed no¬ 
tation suppressing the tensor indices so that ?/„ denotes 
indices at the spacetime point y, etc. 


Following the standard procedure, we define the mul¬ 
tipole moments for the current (106): 


jV^ -y- := (-1)” • • • ay^T’d^^r, (118) 

S(s) 

iVl-VnVO f ...^yr^gvo^^yoy^x' 

S(s) 

(119) 


Then, by applying the lemma (117), we integrate the 
balance equation (107). The result reads 


^ -ai-.-an _ (ai -ao-.-an) i ^ 

—J ” = — nv^ j ' 

ds 

+ V(F£c^d6") - {CcF). (120) 


By boldface symbols we denote the following differential- 
algebraic operators 


^ai...a„b d _ \ ' ^ ai 

^ - 2^ fc! 9 




n + kb d y ... V 


( 121 ) 


OO 1 

= E ^ (122) 


k=0 


which are defined in terms of the multipole moments (C4) 
and (C5). 

The current (106) has a nontrivial structure in that 
it is built from the Noether currents related to the dif- 
feomorphism symmetry of the gravitational theory and 
of the auxiliary objects: vector fields ( and the coupling 
function F. Taking this into account, we can use the 
lemma (117) once again to recast the moments into 


■ai...a„ao ^ ^ai...a„hao ^ao 


(124) 


Here we introduced new operators 

00 - 

k^O 
00 . 

= E ^ • •Va„+,,(126) 

k^O 
00 - 

P^a,...a^bc = ^ 1 y,...a^+ubc^_ .V,„^^,( 127 ) 
fc=0 

defined in terms of the multipole moments of the hy¬ 
permomentum and the energy-momentum currents (Cl)- 
(C3). 

Let us now analyze the moments equation (120). This 
is an infinite system of ordinary differential equations, 
which is common for the equations of motion of extended 
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test bodies. However, one equation is a special one in this 
system. It corresponds to n = 0 and reads explicitly 

^ = Q\^{F£cNdb^) - S{CcF). (128) 

We immediately observe its similarity to (107). More¬ 
over, the structure 

j = P\F(:b) - H\{FVbC) (129) 

is obviously induced by the structure of the current (106). 

After these general derivations, we are now in a posi¬ 
tion to specialize to the pole-dipole case. In the dipole 
approximation, only the first terms proportional to the 
multipole moments up to the dipole order - i.e. 
hab^ j^ab^ J^abc ^ ^abc ^ ^ taken into aC- 

count in (120)-(I29). We then find for the equation (128) 
and for the moment (129): 

^ = q’^/ FC^Nab^ -^JOcF- C^a^cF, (130) 

j = FCa + p'’“V6(ACa) - h\ FVbC. (131) 

Making use of (110)-(112), we immediately see that these 
equations reproduce (115) and (116), respectively. 

VI. CONCLUSIONS 

In this paper we have demonstrated that conserved 
currents can be naturally associated with spacetime dif- 
feomorphisms (represented by vector fields on the space- 
time manifold). Our findings extend previous results 
from Einstein’s gravity with minimal coupling to the gen¬ 
eralized metric-affine gravity theory with a possible non- 
minimal coupling. 

Mathematically, an important role in this construction 
is played by the generalized Killing vector helds, which 
represent the symmetries on the metric-affine spacetime. 
The corresponding formalism is based on the fundamen¬ 
tal geometrical notions of the transposed connection and 
of the Lie derivatives which we heavily used in our pre¬ 
vious studies [27-29]. 

There are many physically interesting applications of 
the results obtained. Among them is the possibility of 
the rigorous computation of the total mass and angular 
momentum for the exact solutions in the gravitational 
theories with and without torsion and nonmetricity. An¬ 
other application is to use the conserved quantity (116) 
to simplify the study of the dynamics of extended mi- 
crostructured test bodies in the generalized gravitational 
field models. Both issues are important for gravitational 
experiments, including space missions, which aim for an 
advanced probe of the geometrical structure of spacetime. 

The results of the present work should be used in the 
context of multipolar approximation schemes, which were 
recently worked out in [68, 79] for a very large class of 
gravitational theories, to further study the dynamics of 
test bodies. Of particular importance is the analysis of 


the detectability of post-Riemannian properties of space- 
time, i.e. the torsion and the nonmetricity. It is worth¬ 
while to mention that the recent literature on the Grav¬ 
ity Probe B experiment [80] encompasses the mislead¬ 
ing claim that the GPB result may set limits on tor¬ 
sion [81, 82]. These erroneous statements were corrected 
in [68, 83-86], demonstrating that in minimally coupled 
gravitational theories the post-Riemannian geometry can 
only be detected with the help of microstructured matter. 

One may notice that the conditions of the generalized 
Killing vectors (79) and (80) imposed on the spacetime 
geometry may be quite strong. This is similar to the 
situation for the manifolds that admit various types of 
collineations. However, even in absence of exact symme¬ 
tries one can consider approximate spacetime symmetries 
along the lines of [56, 87]. We leave the corresponding 
analysis for the future. 
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Appendix A: Notations and conventions 

In order to be consistent with our previous publica¬ 
tions, we choose our main notations and conventions as 
those of [37]. Most importantly, we stick to the defini¬ 
tions of [37] for all the basic geometrical quantities such 
as the curvature, torsion, and nonmetricity, and we use 
the Latin alphabet to label the spacetime coordinate in¬ 
dices. 

The spacetime is modeled as a four-dimensional 
smooth manifold, and its metric has the signature 
(-1-,—,—,—). It should be noted though that our def¬ 
inition of the metrical energy-momentum tensor differs 
by a sign from the definition used in [88]. 

In this work we are widely using tensor densities; we 
denote the densities by the “Fraktur” font to distinguish 
them from the tensor objects. 

Table I displays the list of symbols used in the current 
paper. 
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TABLE I. Directory of symbols. 


Symbol 

Explanation 

Geometrical quantities 

Qab 

Metric 


Determinant of the metric 

a;“, s 

Coordinates, proper time 

c 

Killing vector field 

r.6", Fab^ 

Connection, transposed conn. 

Nab’^ 

Distortion 

Q abc 

Nonmetricity 

rri C 

J-ab 

Torsion 

p d 

-^abc 

Curvature 


Ricci tensor, scalar 

{a^B^ 

Generators coord, transf. 

SS 

Kronecker symbol 

2J 

Gravitational Lagrangian 

^ki. 

Gravitational hypermomentum 


Gener. grav. energy-momentum 


World function, parallel propagator 

Matter quantities 

rj^ab 

Symmetric energy-momentum tensor 


General matter field 

■Cmat 

Matter Lagrangian 

3“ 

Generalized current 


Superpotential 

‘Ik^ 

Canonical energy-momentum 

6‘/ 

Canonical hypermomentum 

n/ 

Spin density 

pf 

Metrical energy-momentum 

Auxiliary quantities 

e 

Infinitesimal parameter 

e' 

Arbitrary vector field 


Multiplet of fields 

I, £ 

General action, Lagrangian 

F 

Coupling function 


Gener. gravitational field momenta 

p -, k -, h -, g ■, 2", 

Test body integrated / generalized 

^ab^ gab^ jab ^ 

moments, velocity 

n -, w“, $ ■, Vt' -' 

Auxiliary variables 

Operators 


Lie derivative 

5 

Transformation under diffeomorph. 

5(s) 

Substantial variation 

di, V^ 

Partial, covariant derivative 

* 

Vi 

Modified cov. derivative 

Vi, Vi 

Gov. density derivative, Riemannian 

44 ?5 

Riemannian quantity 

44 ?? 

Transposed quantity 

“St -,© ■,...” 

Densities “Fraktur” 

H , H , K , 

Auxiliary diff. operators 

p Q 



Appendix B: DifFeomorphism invariance 


The explicit structure of the functions with 

n = 0,1, 2, 3, is as follows 


— di 




~dk9n 


d£ 

dditjj^ 


- Sl£ 


ddiQrnn 

5Z ^ 62, ^ 

dATirT, (Bl) 

OOiTln^ 


dVi 

I r, , A \ i ,B 

M^-512 


+ TTT^- Ok g-n 


6^6)ign 

+a ( 2 


82 

88jgk 


-Qnk + 


dditp^ 

82 , A 


(A \ i I B 
[O- B)k 


88^11}^ 


8Tu3 8r,8 8Tii 


k'^8 


8krin"^ + ^ 8r,rkl^ 


VLk^^ = 


88^21^' 

82 

d8nTi2-' 

A82 

~9kr> 


■ 8nTlk^ - 
282 


88nru"^ 

82 

8dn2lm^ 


8rj:i3n\ (B2) 


88 , 


\i9j)r 


282 

2a£ 




82 


88(^i'ip^ 

282 

m 


-p m 

i kl 


282 


88(^ir\in\'^ 

82 82 


T,J 


5i9*r(j„)'= 88jT(^ni)^ 


(B3) 

(B4) 


Appendix C: Multipole moments 


From the energy-momentum tensor density and the 
hypermomentum density, the integrated multipole mo- 
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merits of arbitrary order, n = 0,1, 2,..., are defined by The integrals are taken over a cross-section S(s) of the 

body’s world tube. Here we introduced 

(C6) 

^v^...v„vov' (C7) 


r,yi---Vn.yo - 


j^y2---yn+iyoyi 


l^y2---yn+iyoyi 


qy3...yn+2yoyiy2 


cyi---yn — 


(- 1 )” f (ci) 

S(s) 

(_ 1 )- f ^y.■■■yr.+^yoy^^^^^ ^ 

S(s) 

(C2) 

(_ 1 )- f ^y2-y.+iyoyi^^^^ ^ 

E(s) 

(C3) 

(_ 1 )- f ^,y3-yr.+.yoyi^^^^gy2^^ ^ 

S{s) 

x6’"'>"’i"="w"’3dS^3, (C4) 


E(s) 


(C5) 


In the derivation of the equations of motion we made 
use of the bitensor formalism; see, e.g., [78, 89, 90] for 
introductions and references. In particular, the world 

/ y 

function is defined as an integral a{x,y) := I f ds 

\x 

over the geodesic curve connecting the spacetime points 
X and y, where the upper/lower sign is chosen for time¬ 
like/spacelike curves, respectively. Note that our curva¬ 
ture conventions differ from those in [78, 90]. Indices 
attached to the world function always denote covariant 
derivatives, at the given point, i.e. Uy := Vycr; hence, we 
do not make explicit use of the semicolon in the case of 
the world function. The parallel propagator by gyx{x, y) 
allows for the parallel transportation of objects along the 
unique geodesic that links the points x and y. For exam¬ 
ple, given a vector at x, the corresponding vector at y 
is obtained by means of the parallel transport along the 
geodesic curve as = gyx{x,y)V^. For more details 
see, e.g., section 5 in [90]. A compact summary of useful 
formulas in the context of the bitensor formalism, as well 
as a review of the multipolar formalism employed here 
can also be found in [79]. 
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